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Abstract 

Spacecraft control, state estimation, ard fault -detection-and- 
isolation systems are affected by unknown ' ariations in the vehicle 
mass properties. It is often difficult to accurately measure inertia 
terms on the ground, and mass properties (an change on-orbit as 
fuel is expended, the configuration changes or payloads are added 
or removed. Recursive least squares -based algorithms that use gyro 
signals to identify the center of mass and inverse inertia matrix are 
presented. They are applied in simulation to 3 thruster-controlled 
vehicles: the X-38 and Mini-AERCam under development at 
NASA-JSC, and the S4, an air-bearing spacecraft simulator at the 
NASA-Ames Smart Systems Research Lab I SSRL). 

1. Introduction 

The mass-property identification (ID) algorithms presented here 
were developed through application to two thruster-controlled 
spacecraft presently under development at NASA Johnson Space 
Center: the X-38 [14] and the Mini-AERCam, and to the SSRL S4. 


♦ 



Fig 1: X-38, Mini-AERCam, and SSRL air-bearing vehicle, S4 


The goal is to identify, on-line, the mass property parameters of a 
thruster-controlled spacecraft. Specifically, he terms of the inverse 
inertia matrix and the location of the mass center are identified 
using gyro measurements during periods of thruster firing. The 
application vehicles, X-38, Mini-AERCam, and S4 are all 
exclusively thruster controlled and the deri vation follows that, but 
it is a simple extension to accommodate CMGs or reaction wheels 
(as in [2]). Accelerometers would improve the ID, but gyros only 
are used, broadening the applicability of the resulting algorithm. 
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1.1 Related research 

Tanygin and Williams [9] developed a least squares (LS) based 
algorithm to identify mass properties for a spinning vehicle during 
coasting maneuvers. Bergmann et al [2][3][8] developed a Kalman- 
filter based approach to the mass-property ID problem. A careful 
comparison of this approach to the recursive least squares (RLS) 
approach presented here has not yet been completed. Wilson and 
Rock [1 1 ] [ 1 2] developed an ID method based on exponentially 
weighted RLS using accelerometer and angular rate sensors. The 
acceleration created by each thruster (reflecting both mass and 
thruster properties) was identified. A neural network then provide 
adaptive control reconfiguration to multiple destabilizing hard and 
soft thruster failures. This was implemented on a ^degree of- 
freedom air-bearing vehicle. Wilson, et al developed a model-based 
thruster fault detection and isolation (FDI) system for the vehicles 
presented here [13]. These mass-property ID algorithms were 
developed to improve model accuracy for that system. 

1.2 Least-squares identification 

In this paper, vehicle mass properties are identified using LS 
methods in which the sensor data is fit to the underlying equations 
of motion (EOM) such that the identified parameter values 
minimize the squared erTor (where error is, for example, sensor 
data minus the ideal sensor data that would occur with zero noise 
and using the identified parameter values). 

The standard form for a linear least squares problem is given as 

Ax = b + £ or Ax = b 

where b is a vector of noise-free measurements, E is a vector of 
measurement noise, X contains the parameters to be identified, 
and matrix A contains known variables and parameters (i.e. , A is 
noise-free). The = in Ax = b indicates that the left and nght 
sides of the equation would be equal if noise were not present [6]. 

The LS ID solution, X , minimizes the sum of the squares of the 
error, Ax — b . If the problem at hand can be put into this standard 
form, X can be solved directly using a batch algorithm, 

x.L'Af'Sb. or an exactly equivalent recursive algorithm 
[4][5][6]. Manipulating the original equations into the form 
Ax = b so that the standard LS solution can be solved is often 
the primary challenge, requiring careful, application-dependent 
decisions regarding approximations. 

The rotational EOM (derived later) contain all of the parameters to 
be identified: center of mass (CM) location, contained in L , and 
inertia. Unfortunately, these parameters multiply one another, and 
cannot all simultaneously be manipulated into the desired linear 

form, Ax = b . The approach taken here addresses this issue by 
segmenting the ID problem into two sub-problems which both 
allow closed form solution. 


2. Approach 

For the vehicle applications studied, the presence of disturbance 
torques, imperfect rate measurements, and significant pulse-to- 
pulse thruster variability biased development towards one of 
minimal complexity (e.g. f dropping higher order terms, segmenting 
the problem, etc.). Before going into the equations, it is useful to 
understand the basic physics behind the ID approach: 

1. The rotational EOM mathematically describe how forces and 
torques will affect rotational accelerations. The inertia matrix 
governs the acceleration that will result from an applied net 

torque, allowing LS ID of I~ l . The CM location determines 
the angular acceleration that will result "rom an applied force. 

2. When thrusters are fired to produce a pure torque only, the 
CM location has no effect on angular acceleration. So the CM 
ID cannot be updated using this data. Also, inertia ID will not 
be biased by an incorrect CM ID estimate. 

3. When force and torque are applied simultaneously, inertia and 
CM can each be identified individually. 

4. When no external forces or torques an applied, disturbances 
due to thruster variability are removed. With sufficiently high 
angular rates and sufficiently low sens< >r noise, accelerometer 
signals can be used along with the gyros to ID the CM 
location. This has been developed and implemented by the 
authors for a vehicle with relatively lew sensor noise levels 

[9]), and is not presented here. 

Specifics regarding this approach for gyro-bi sed mass-property ID: 

• Due to the form of the EOM, the invene of the inertia matrix 

is identified (i.e., I 1 vs. 1 ). 

• Only the six independent elements of I 1 are identified. 

• / _1 and CM updates are made only when thrusters are firing. 

• Perfect knowledge of thruster failures and biases is assumed. 

• If pure torques occur frequently, item 2 above is followed, so 

/ -l updates are made only when a pure torque is present 
(assuming that fewer clean data points are better than more 
noisy ones, and that the CM error would bias the estimate). 

3. Derivation of mass-property RLS ID algorithms 

In this section, the relevant EOM are derived, containing the 
measurements and parameters of interest. Then these EOM are 
manipulated into forms that allow ID, first of the mass center and 
second of the inverse inertia matrix. Once the equations are in the 
proper form, either batch LS or RLS can be implemented. 
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3.1 Equations of motion 

Starting with Euler’s dynamical equatior , and assuming the 
spacecraft inertia matrix is constant, the rotat onal EOM are [1]: 
c q - /-i( T -(olco) > or, showing individual elements, 
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where I is the spacecraft inertia matrix (and I- , indicates the 

(i,j) element of / ' rather than 1 / 1 tJ ), CO is the angular 

velocity of the body -fixed frame with respect to an inertial 
reference frame, T is the sum of all torques on the body (thrusters, 


reaction wheels, CMGs, disturbances, etc.), and (0 is used to 
represent the matrix-multiply implementation of the cross product. 

The total torque, T , about the true CM (which is not known 
exactly) due to the thrusters and torque disturbance is modeled as: 

T = T thruster A ^ disturb = (^ X D ) F t ^disturb 

where T lhrusters , represents the net thruster torque. L and D are 

[5-by- ] matrices (N is the number of thrusters) containing the 
location and direction of each thruster, measured in the body frame. 
Each column contains the information for an individual thruster, 

and the cross product operator shown here indicates that the i 

column of L> is crossed with the / column of D and entered 

into the i th column of the resulting [J-by- N ] matrix. 

F k , the thrust magnitude from each individual thruster at each 

control update, k , can be expressed as the following function, 
which accounts for (in order in the equation): reduction in thrust 
due to blowdown; nominal thruster magnitude; constant random 
bias added to the nominal value; random pulsed o-pu l se offset 
added to the nominal value; thruster firing (on/off). A further 

refinement in calculating F k would be to also model the transient 
nature of the thrusters, accounting for rise, trail off, and latency. 

F t =B(F nom + F bios +F rM )T k 

B and F bias can be identified separately. In this analysis, the ID 
algorithm uses their nominal values (B= 1 ,F bias =0) * but their true 
random values are used in the simulation for testing. 

Combining the above results, the full rotational EOM reduce to 

ib=r l ((LxD)B(F om + F b , a , + F^)T k +T^ urb -a>x(fa») 

The EOM are modeled assuming that the acceleration is constant 
during each control sample period; the ID algorithms require 
acceleration estimates corresponding to these control sample 
periods. The applications presented here estimate angular 

acceleration, (O , using rate gyros [13], but any method could be 
used (e.g., using Kalman filtering, star trackers, accelerometers). 

3.2 Identification of the mass center 

The mass center, C (measured in the structural frame), determines 
the origin of the body frame, and thereby determines the value of 

L , which contains the locations of each thruster in the body frame. 
Similarly, A , the difference between actual and nominal values of 
C determines L . A is the value that will be identified here. 
CaC„,+4;i=i.„-A(l 1 ••• l] 

Recalling the rotational EOM, substituting and rearranging some 
terms, and setting B=-\,F bmi =0, F randmk = 0,f Jis , ur6 = 0 , 

(b=r'{(LxD)F HO J t -WY.(Ia))) 

Id) +fflx(to) = - A[l 1 ... 1 ]>xD)F nom T k 

Trying to get this equation in the form Ax = h with A as X , 

Id) + cox (la)) = (L nom x D)F no J' k - (A[l 1 ... l]xD)F nom T k 


Since the columns of A[l 1 ... l] are the second 

term on the right can be regrouped. Note as a counter-example that 
the first term on the right cannot be similarly regrouped because the 


columns of L nom are not all the same, 


Id) + o)x(Io)) =(L Mm x D)F nom T k - ^x(DF no T k )) 

Ax (DF mem T k ) = (L nom x D)F mm T k - Id) -O)x(Iw) 

By the anti-commutative property of the cross product (that 
a xb = -(bxa))) the left side is changed to 

(DF no T k )x A = (L nom x D)F ro T t - !ti) -0) x (to) 


The right side (the b in Ax = b ) can be calculated. Introduce a 
variable, c k = DF nom T k (a 3-by-l vector), then re-write using the 
matrix- multiply implementation of the cros^ -product, as 
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= (L nom x D)F a J k - Id) -o)x(h y) 


The measurements in this equation show ip in the (D and (D 
terms only. The (D X {1(D) term is very srr all (and can be omitted 

for the applications studied here), so CO is the main noise 
contributor. The LS problem is defined (and closed-form solution 
developed) assuming the measurement noise enters without pre- 
multiplication. So, p re-multiplying by 1 1 and changing sign, 
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(b + I \co x(I(o)) - / 1 {L nom x D)F nor T k 
The a) X(I(0) term is still pre-multiplied, but that is likely to be an 

insignificant number. So at each upda e, c k = DF nom T k is 
calculated and plugged in to the above re-formulation of the 
equation of motion, which conforms to Ax = b , where 
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b k = (i) + F l (co x (/ c o )) - r 1 (L nom x D )F nom T k 


3.3 Identification of the inverse inertU matrix 

ID of I' 1 is similar to that for A: the rotational EOM are 
manipulated to a form as close as possible 'o Ax = b . Recalling 
the original rotational EOM and simplifying as before, 

(i) = r l ((LxD)F no T k -cox(Ia))) 

I " l ((L x D)F no T k -wx( Ico)) = d) 

Introduce a variable, a k = (L x D)F nom T k -cox (Ico) > so 
l~ x a — (i ) . The /“' a k matrix-vector muliiplication is converted 
into an equivalent expression where the independent terms in 
/“* appear in a vector that is pre-multipliec by a matnx. 
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So now the equation is in the standard form, Ax — b , where 


Au = 


a 2 a 3 
a \ 




In 
1 22 
In 
In 
In 


b k =(b k 


For both mass-center ID and inertia ID, measurements at multiple 
time steps are combined as follows if the batch LS solution is used. 
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3.4 Deviations from standard LS ID form 

The recursive or batch LS solution, x = (a t A) A 1 b > will truly 
minimize the quadratic error if the manipulated equations comply 
with the standard form, Ax=b + £ . Deviations from this form in 
the results of the preceding manipulations are that: 

1 . Noisy measurements appear in the (0 X ( 1(D) term in the A 
matrix. However, for the relatively slow rotational speeds in 
many spacecraft applications, this term is negligible. 

2. Other terms in A and b are not known perfectly: L , D , 
B F. etc. are all estimated or nominal values. 

* bias ’ 

3. Random variables F mndom and X disturh that were set to zero 
do not appear directly in the £ term as they should. 

4. CM ID uses nominal or estimated values for / and 1 
Inertia ID uses nominal or estimated values for L (CM). 

Due to the form of the underlying EOM, it is not possible to exactly 
comply with the standard form. The preceding manipulations and 
approximations to the EOM were chosen for this class of 
applications, attempting to minimize the expected effects of the 
inevitable resulting deviations. With different typical values of (D , 
disturbances, sensor accuracy, actuator variability, control policy, 
etc. y different formulations may be better (e g., [9]). 

Also, the assumptions that acceleration is constant during each 
control update, and that thruster response time is zero will degrade 
the ID results if not compensated for. 


4. Simulation Results 

A dynamic simulation was developed using MATLAB [7] that 
included: the rotational and translational EOM; different sample 
rates for control and sensing; a thruster control system; and realistic 
variations in all relevant system parameters (mass properties, 
thruster properties and pulse-to-pulse noise, sensor noise). 

The accuracy of the ID results depends on the sensor noise, thruster 
variability, and variability in non-identified system parameters 

(such as thruster direction, D , and bias, F pias )• These RLS mass- 
property ID algorithms have been implemented and successfully 
tested for the three dissimilar spacecraft mentioned. In this section, 
results for one, the Mini-AERCam, are presented as an example. 

The Mini-AERCam (upper right of Figure 1) has 12 thrusters 
controlling all 6 degrees of freedom, and uses MEMS gyros. In this 
test, the system was excited automatically by randomly perturbing 
the desired position and attitude, causing thrusters to fire, providing 
both torques and net forces on the vehicle. 



Fig 2: Mini-AERCam inertia-inverse ID results 

Figure 2 shows the results for ID of the inverse inertia matrix. True 
values are drawn as solid lines; ID updates a^e drawn as dots. Since 
pure torques are common for this vehicle, irertia ID updates occur 
only when nominally pure torques are applied - resulting in the 
gaps between updates. 
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Fig 3: Mini-AERCam mass-center ID results 

CM ID results are shown in Figure 3. As with inertia ID, the ID is 
initialized at the nominal values, which is zero in this case. CM 
information is not present when pure torques are applied, so there 
are no updates at these times (or when no thrusters fire), resulting 
in the gap between 2 and 5 seconds. In both inertia and CM ID, ID 
accuracy is better than the vehicle spec based on ground analysis 
and measurement (±5 mm for CM in this example), making the ID 
useful for on-line adaptive control. 

Although they are run independently in thej e tests, use of CM ID 
results in the inertia ID calculations (i.£., losing ID ed instead of 
nominal CM information) and conversely wc uld improve results. 


5. Conclusions 

Separate mass-center and inverse- inertia-matrix ID algorithms have 
been developed and applied in simulation to the X-38, Mini- 
AERCam, and S4 thruster-controlled vehicles. Using gyro signals 
only, and based on recursive least squares, the algorithms reliably 
and accurately ID mass properties for these vehicles in the presence 
of several significant noise sources. The algorithms are 
computationally efficient and can be run either on-line for adaptive 
control (RLS) or off-line for post-flight analysis (batch LS). 
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